Gossamer Superconductivity near Antiferromagnetic Mott Insulator in Layered 

Organic Conductors 
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Layered organic superconductors are on the verge of the Mott insulator. We use Gutzwiller 
variational method to study a Hubbard model including a spin exchange coupling term. The ground 
state is found to be a Gossamer superconductor at small on-site Coulomb repulsion U and an 
antiferromagnetic Mott insulator at large U , separated by a first order phase transition. Our theory 
is qualitatively consistent with major experiments reported in organic superconductors. 
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There has been much interest recently on the novel 
physics of layered organic superconductor 0, S US . 
These compounds share most common physical proper- 
ties with the high-T c superconductor but typically with 
much reduced temperature and energy scales. k-(BEDT- 
TTF)2X (X=anion) is a family of the best characterized 
organic superconductors, where the quasi-2-dimcnsional 
(2D) Fermi surface has been observed and a direct first 
order transition between antiferromagnetic (AFM) in- 
sulator and superconductor can be tuned by applied 
pressure or magnetic fields 0, S S IS E3- The re- 
semblance of its pressure-temperature phase diagram to 
that of the carrier-density-temperature phase diagram 
in cuprates and the fact of close proximity between the 
superconducting (SC) and AFM insulating phases have 
been taken as evidences for similar mechanisms govern- 
ing high T c superconductors. There have been strong evi- 
dences that the organic superconductors are at the verge 
of the Mott insulator [(| 0, exhibiting the pseudogap 
phenomenon @ . While an ongoing debate persists as to 
the precise symmetry of the singlet pairin g, m ore recent 
NMR [U [?|, angular dependent STM [13| and ther- 
mal conductivity measurements in the vortex state EH 
indicate a d x 2_ y 2 symmetry. 

The low energy electronic structure of the organic su- 
perconductors is well approximated by a 2D Hubbard 
model at the half filling. Different from the cuprates, 
organic compounds can be SC at the half filling, which 
makes t-J model inappropriate to describe its SC state. 
Most theoretical works so far have taken a weak-coupling 
approach, in which a Hartree-Fock mean field | 15l| . 
fluctuation-exchange approximation IT^. ITtI ITsl Il9| or 
random phase approximation method |20| are used. The 
weak coupling theory gives a phase diagram of the AFM 
and SC states qualitatively consistent with the experi- 
ments. However, the weak coupling theory has difficul- 
ties to address the Mott insulator or the pseudogap phe- 
nomenon [ESQ- The transition between SC and AFM 
has also been investigated by using renormalization group 
method [2l| . 

Very recently, Laughlin has proposed a Gossamer 



Hamiltonian of which a partially Gutzwiller projected 
BCS state is an exact ground state with a tiny super- 
fluid density at the half filling [z| . In that Hamiltonian, 
the SC state has an instability toward the AFM order- 
ing [2i| . Some of the present authors have exam- 
ined the Gossamer superconductor, the Mott insulator, 
and the resonating valence bond (RVB) state pr j |27[ Esf 
in strongly correlated electron systems with the hope to 
unify the superconductivity in cuprates and in organic 
compounds [29j. In our previous study, we focused on 
the metallic/SC and insulating nature of the problem and 
neglected the antiferromagnetism. A related approach 
was recently taken by Baskaran |30l |. who introduced a 
two-species t-J model to describe independent motions 
of empty sites and doubly occupied sites in an otherwise 
spin-1/2 background, and discussed the relevance of the 
model to the organic superconductors. 

In this Letter, we use Gutzwiller's variational method 
to study the interplay between SC and AFM states in a 
modified Hubbard model in 2D given by Eqn. (1) be- 
low. By using a renormalized mean field theory devel- 
oped early for the t-J model [27|], we find that at the 
half filling the ground state is an AFM Mott insulator 
at large on-site repulsion U and a Gossamer supercon- 
ductor at small U, followed by a normal metallic state 
at further smaller U. The transition between the AFM 
and SC phases is first order, and there is no co-existence 
of the two phases at the half filling. The doping depen- 
dence of the model at large U is similar to that of the t-J 
model 26] . Our results are qualitatively consistent with 
major experiments in organic superconductors. 

We consider a modified Hubbard model on a square 
lattice, 

H = U^n^riii - *«( c L c jV + h - c ) 

i (ij)o 

+j §i ■ §j — n nig (1) 

(ij) iff 

In the above Hamiltonian, Ci„ is an annihilation opera- 
tor of an electron at site i with spin c, rii a = cj^cv, and 
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U > is the on-site Coulomb repulsion. The non-zero 
hopping integrals are fcy = t for the nearest neighbor 
(n.n.) pairs and £y = t' for the next n.n. pairs along 
[1,1] direction. Si is a spin-1/2 operator, and the summa- 
tion in the spin exchange term is over all the n.n. pairs. 
We consider 1 > t' /t > 0, suitable for the organic com- 
pounds. This Hamiltonian contains an additional spin 
exchange term to the standard Hubbard model. In the 
limit U ^> t, the model is reduced to the t-t'-J model. At 
the half filling, the large U limit of the model is reduced 
to the AFM Heisenberg model with an AFM ground state 
at small values of t' /t. At small U, we expect a metallic or 
a SC ground state. We believe that the model combined 
with the Gutzwiller trial wavefunction approach, Eqns. 
(2-3) below is appropriate to study the phase transitions 
in organic superconductors. Note that the direct applica- 
tion of the Gutzwiller trial wavefunction to the Hubbard 
model is hardly to obtain the SC pairing because of the 
non-explicit form of the spin-spin exchange interaction in 
the Hamiltonian. 

To study the phase transition between the AFM and 
SC states, we consider a partially Gutzwiller projected 
spin density wave (SDW)-BCS wavefunction |3ll l32|. 

I*gs> =n( 1 - Q ^T^l)l*o) (2) 

i 

\y°) = U( u k+ v A ] d U i )l0) (3) 

where d %a = cos{ 9 -f)c ia - asm {^)c^ a , and Q = 
(-7T, 7r) is the magnetic wave vector. J}(1 — an^n^) is a 

i 

Gutzwiller projection operator, which partially projects 
out the doubly occupied electron states on every lattice 
site and < a < 1 measures the strength of the pro- 
jection. Obviously, a = 0, and a = 1 correspond to 
a non-projected and a completely projected states, re- 
spectively. At Oj: = 0, we have dz = cj ct , and | ^>gs) 
is reduced to a partially projected BCS state, which we 
shall loosely call it Gossamer SC state [2^, ■ In the 
limit uj:vj: = 0, | ^o) is reduced to a SDW state. The 
variational parameters arc «£, v^, 9^ and a. Such a wave- 
function should enable us to study the phase transition 
between the AFM and SC states. The metallic or insu- 
lating phase can be determined by the continuity of the 
chemical potential. 

To carry out the variation, we apply the Gutzwiller 
approximation to replace the effect of the projection op- 
erator by a set of renormalization factors, which are de- 
termined by statistical countings [53, l33l lijil l3q . Let 
(O) be the expectation value of the operator O in the 
state I^gs), and (O)o be that in the state l^o)- The 
Gutzwiller approximation gives 

< C L C ^> = 9t (4t c j CT )o, {S t ■ Sj) = g B {Si ■ Sj) (4) 

where g's are determined by the ratio of the probability 
of the corresponding physical processes in the projected 



and unprojected states [23 . We introduce a sublattice 
magnetization for sublattices A and B, 

mo = \i n A\ - n A i)a = -^(nst _ n Bi)o (5) 

g's arc then functions of the electron density n, too, and 
the double occupation number d = (n^n^), 

g s = (n — 2d) 2 / {n — 2n + n_) 2 

g\ 3 = G l G 3 

G A = 3 y 4 [ S (l-n_) + v /Z27^] 

G B = 5 y 4 [ S (l-n + ) + v ^7rfAT] (6) 

In the above equations, n± = || ± too, and s = 

\J (i-n7)"u-» ) • ^ ne superindex in G refers to the sub- 
lattice of the site. Note that there is a one-to-one corre- 
spondence between d and a given by, 

1 — a = s 2 d/g s n + ri- (7) 

In the absence of the sublattice magnetization, g t and g s 
in Eqn, (6) are reduced to their values in the uniform 
state |25j, which are further reduced, in fully projected 
case (a — 1 or d — 0), to the values in the RVB state |27| . 
Within the Gutzwiller approximation, the variation of 
the projected state for H in (1) is reduced to the vari- 
ation of the unprojected state | \&o) for a renormalized 
Hamiltonian H e ff, 

H eff = Ud- gl J (ct a c ja + h.c.) 

+g s J} j Sj-Sj - /J,}^n i(T (8) 

(ij) ia 

To proceed further, we introduce a self-energy \ and a 
d-wave pairing amplitude A, 

x = ^2(4^+^)0 = J2(4 a c i+ya ) (9) 

(X (T 

A = y^(crc. ta -Ci +a _ er )o = - y^^(<TCi a c t+y - a ) (10) 

(X (J 

The singlet SC order parameter A$c ~ 9scA, with 
9sc = (gf A + 9t B )l^- The pairing amplitude and 
the SDW state described below defines the variation of 
| ^o)- As in the usual SDW variation, we choose cos#£ = 

e fe/£fc' where e fe = ~( 2t 9f B + 3J5 s x/4)7fe !+ is the ki- 
netic energy including a self-energy term of x, and £g = 

yJel + A* AF (k), with A AF (k) = A af + t'(g AA - gf B )C £ . 
A a f is a variational parameter to determine to . The sec- 
ond term in Aaf arises from a spin-dependent hopping 
process along the [1,1] direction in H e ff. In the above 
equations, we have denoted 7^ ± = cosk x ± cosk yi and 
= cos (k x + k y ). With the above variational wavefunc- 
tion, we calculate the expectation value of H e f f and find 
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FIG. 1: Pairing amplitude A, SC order parameter Asc, and 
AFM order parameter m (top panel), and electron double 
occupancy number d and the projection parameter a (bottom 
panel), as functions of U for J/t = 0.5 and t /t = 0.8. 



the ground state energy. 

E = Ud-Ag t t X + g? A (H+) a 

~(3g s J/4)(A 2 + X 2 )-2Jg s m 2 



(11) 



where too, Xi n i an d A are the solutions of their corre- 
sponding self-consistent equations. The two additional 
variational parameters d and A a f are to minimalize the 
ground state energy. Note that < d < d , with 
d = (n^ni^Q. In Eqn. (11), (Hp) are given by 



H ?)o = -^ECfe[^(lTsin2%) 



,2 +Q -(l±sin2^)] 



where the summation of k runs over the reduced Brillouin 
zone, and 



i 



W fc+Q 



1 



1 + fe 



with _5T = 
fc 



A|, and A s 



(3/4)J 5a A 7S) _, 

M = ^ + <V A + .Oc fe -. 

We are now ready to discuss our results. We shall 
mainly discuss the half filled case. At the half filling, 
there is a critical U c to separate a metallic or SC state 
at a small U from an AFM insulator at a large U, and 
the transition is first order with no co-existence of the 
two phases. These features are demonstrated in Fig. 1. 
There are two regimes in U. At U < U c (^ 5.5 1), to = 
while A and Asc increase monotonically as U increases. 
Asc is slightly smaller than A. This is a SC state with- 
out AFM ordering. At U > U c , A = A S c = 0, while 
to = ^/g s mo changes abruptly from zero at U < U c to a 
saturated value of 0.45. We have calculated the chemical 
potential around the half filling and found it is discon- 
tinuous in the AFM state so that it is an insulator. As 
we can see from the bottom panel of Fig. 1, as U in- 
creases, d decreases with a sudden drop at U c indicating 
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FIG. 2: Left panel: phase diagram of t' v.s. U for J/t — 0.5 at 
the half filling. The arrow indicates the flow of the parameters 
under the pressure. Right panel: schematic phase diagram of 
organic superconductors. 



the electron's localization in the insulating phase, and 
a increases to its maximum in the SC phase followed 
by a discontinuous drop to zero at U — U c . The lat- 
ter indicates the absence of the projection in the AFM 
phase so that we have m = m .'Hi . We have also cal- 
culated these quantities with different values of J/t and 
t' /t and the results are qualitatively similar except that 
A becomes very tiny at smaller J/t. Our results are con- 
sistent with major experiments in organic superconduc- 
tors. As shown in the pressure experiments pj, the phase 
transition is first order and the phase boundary between 
the AFM and SC states merges with the phase boundary 
between the insulating and metallic states. In our the- 
ory, the AFM state is always a Mott insulator. Recent 
NMR experiments 0, |1] show the proximity of pseudo- 
gapped superconductor and a commensurate AFM or- 
dering with a finite moment of OA/is (or 0.26/^s) for k- 
(BEDT-TTF) 2 Cu[N(CN) 2 ]C1 or Br at low temperatures, 
which suggests that the magnetic ordering is driven by 
electron's strong correlation rather than by the Fermi 
surface nesting. In the Gossamer SC state, the quasi- 
particle energy is governed by A [2{| |2(| |23, which is 
larger than the SC order parameter, implying a pseudo- 
gap phase [3]]]. The small difference between Asc an d 
A in our theory is partly due to the phenomenological 
model we use, which more favors AFM state than the 
Hubbard model does at moderate or large U. We expect 
the phase boundary in a more accurate theory will be 
shifted to the larger U and Asc/ A will be smaller. 

Fig. 2 displays the phase diagram in the parameter 
space of t' and U with fixed J/t — 0.5 at the half fill- 
ing. There are three distinct phases. The system is in 
the AFM phase at large U and small t', the paramag- 
netic metallic phase at small U and large t', and the SC 
phase at the intermediate parameter region. Here we 
have defined a paramagnetic metallic phase if A < 0.01. 
At this very small A, the energy difference between a 
SC state and a normal metallic state is practically in- 
distinguishable. The phase boundary between the SC 
and normal states thus obtained is indicated by a dashed 
line |38j . For comparison, a schematic phase diagram 
abstracted from experimental measurements is shown 
at the right panel. Details of the pressure-temperature 
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FIG. 3: Doping 5 dependence of A, Asc and m (top panel), 
and d and a (bottom panel) for U/t = 8, J/t — 0.5 and 
t'/t = 0.8. 

phase diagram of the AFM insulating salt have been re- 
ported HUGS- 

The effect of pressure in the schematic 
phase diagram is to decrease U/t or to increase t'/t. Our 
theory is consistent with the general features of this ex- 
perimental phase diagram. Note that the t-U-J model 
does not represent the Hubbard model at small U. As 
it is well known, the ground state of the Hubbard model 
with t 1 = at the half filling is an antiferromagnet. In 
our study of Eqn. (1), J is considered as an independent 
parameter, so that the J-term together with the kinetic 
term is in favor of a metallic state at U = or small U/t. 

Away from the half filling, similar calculations are con- 
ducted. The doping (5 = 1 — n) dependences of various 



quantities are plotted in Fig. 3 for U > U c . As S in- 
creases from zero to a critical doping S c « 0.07, the AFM 
order parameter m decreases slightly, and Asc increases 
initially, then saturates. The ground state is an unpro- 
jected SDW and SC state (a = 0), where the AFM and 
SC phases co-exist. At (5 = S c , m drops to zero. The sud- 
den disappearance of the AFM order strongly enhances 
the SC pairing. Asc has a jump at S = 5 C followed by a 
slow decrease as 5 further increases. In the region S > 6 C , 
we have a pure d-wave SC state. The essential physics 
here is similar to the doped RVB state j2||, except that 
here we have a strong first order phase transition on the 
AFM ordering at S c , a point which requires further study. 

In summary, we have presented a strong coupling varia- 
tional theory to examine the superconductivity near an- 
tiferromagnetic Mott insulator in layered organic con- 
ductors by using a Hubbard model including a spin-spin 
coupling term. The theory appears qualitatively consis- 
tent with a number of major experiments, such as the 
first order phase transition between AFM Mott insula- 
tor and superconductor under pressure, the merge of the 
metal-insulator transition and the AFM-SC transition 
point, the pseudogapped phenomenon in the SC state, 
the large magnetic moment in the AFM phase, and the 
transition to the normal metallic phase at high pressure. 
The present theory may be further improved by using 
variational Monte Carlo calculations on the expectation 
values and by developing a more accurate Hamiltonian- 
trial-wavefunction approach describing the physics of the 
layered organic conductors. 

This work is partially supported by Chinese Academy 
of Sciences, and by the RGC of Hong Kong. 
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